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INTRODUCTION 


The  time  and  frequency  domains  are  equivalent  displays  of  seismic  trace 
information,  though  some  qualities  of  the  signal  are  more  easily  observed  in 
one  domain  than  the  other.  The  relative  frequency  excitation  of  Lg,  for 
instance,  is  most  easily  viewed  in  the  frequency  domain,  but  such  waveform 
qualities  as  the  sequence  in  which  pulses  arrive  in  the  wave  train  or  the 
sharpness  of  pulse  onset  are  most  easily  studied  in  the  time  domain  (Murphy 
and  Bennett,  1982;  Burdick  and  Helmberger,  1988;  Blandford,  1981).  Because  of 
the  tremendous  complexity  of  high  frequency  regional  data,  most  attempts  at 
using  it  for  discrimination  purposes  have  involved  analysis  of  the  frequency 
content  of  the  various  arrivals  either  through  transforming  selected  windows 
or  through  multiple  bandpass  filtering.  In  this  report,  we  examine  the 
alternative  and  attempt  to  discriminate  events  using  those  waveform 
characteristics  most  easily  observed  in  the  time  domain. 

The  greatest  difficulty  in  interpreting  short  period  regional  data  is 
that  it  is  so  strongly  affected  by  scattering  phenomena.  The  best  approach  to 
analyzing  it  may  be  to  use  array  processing  techniques  on  data  recorded  at 
technically  advanced  seismic  arrays.  Such  techniques  measure  the  average 
properties  of  signals  arriving  in  given  velocity  windows.  In  this  work,  we 
are  attempting  to  characterize  the  average  properties  of  some  of  the  more 
prominent  regional  phases  using  fundamental  array  processing  techniques. 
However,  we  have  reversed  the  role  of  sources  and  stations.  We  use  an  array 
of  sources  recorded  at  a  single  station  rather  than  the  converse. 

In  analyzing  a  test  discrimination  data  set  from  the  western  U.  S.,  we 
have  discovered  that  che  onset  of  Pn  is  always  very  similar  for  explosions  and 
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that  few  earthquakes  have  this  unique  waveform  character.  We  have  also 
discovered  that  the  regional  phase  Pg  is  actually  composed  of  a  sequence  of 
sub-arrivals  which  correspond  to  successively  higher  orders  of  reverberation 
in  the  crust  (Burdicx  and  Helmberger,  1988).  The  scattering  associated  with 
each  new  reverberation  accounts  for  the  long  duration  and  complexity  of  p 
We  suspect  that  an  analogous  set  of  reverberations  with  scattering  may  account 
for  Sg.  In  analyzing  these  crustal  resonance  phases,  we  found  that  the  first 
reverberation  is  dominated  by  PmP  type  rays,  but  the  second  is  not  dominated 
by  second  multiples  of  P  energy  as  might  have  been  expected.  Rather,  it  is 
dominated  by  waves  with  at  least  one  converted  leg.  In  viewing  high  frequency 
data,  we  have  found  that  the  second  resonance  phase  has  markedly  different 
particle  motion  than  the  first  and  that  the  second  phase  is  much  richer  in  S 
energy.  This  suggests  that  discrimination  may  be  possible  through  analysis  of 
the  evolution  of  particle  motion  of  Pg  with  time.  We  first  discuss  the  Pn 
discriminant  and  then  the  Pg. 
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THE  P„  WAVEFORM  DISCRIMINANT 

The  situation  of  having  data  available  from  many  sources  very  close  to 
each  other  which  are  similar  in  character  such  as  NTS  explosions  is  a  unique 
one.  The  timing  and  locations  of  the  events  are  known  exactly.  This  turns 
any  single  regional  station  in  the  western  U.  S.  into  the  equivalent  of  a 
regional  seismic  array.  Any  of  the  standard  array  processing  techniques  can 
be  used  with  the  role  of  sources  and  stations  being  reversed.  There  is  some 
variability  in  source  time  history  and  near  source  structure,  but  on  the  other 
hand  the  receiver  structure  is  constant.  This  fact  is  being  exploited  in  the 
studies  discussed  in  the  following  in  an  attempt  to  find  what  properties  of 
regional  data  are  stable  enough  on  average  to  be  used  for  discrimination 
purposes.  There  are  a  number  of  high  quality  regional  seismic  stations  around 
NTS  from  which  to  choose.  These  include  the  DWWSSN  stations  ALQ  and  JAS ,  the 
four  broad  band  LLNL  stations  and  new  stations  at  Pasadena  and  Pinyon  Flat. 
These  last  two  are  instrumented  with  the  state  of  the  art  Streckeisens  which 
are  installed  at  Grafenberg  and  will  be  installed  at  the  new  international 
stations  in  the  Soviet  Union.  For  this,  the  exploratory  part  of  the 
investigation,  we  simply  use  those  stations  for  which  the  data  are  routinely 
included  on  the  USGS  Network  Day  Tape.  These  are  the  WWSSN  stations  ALQ  and 
JAS. 

To  begin  the  study,  data  bases  of  explosion  and  earthquake  signals  were 
assembled  from  the  day  tape.  The  short  and  long  period  signals  were  plotted 
n  expanded  time  scales  and  examined  for  features  that  tended  to  discriminate 
r’e  two  source  types.  In  so  doing,  it  was  noted  that  the  first  few  swings  of 
P  were  very  consistent  in  shape  and  contained  a  very  interesting  feature. 
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Six  typical  Pn's  from  Yucca  events  as  observed  at  ALQ  and  six  others  at  JAS 
are  shown  in  Figures  1  and  2.  The  variation  in  event  size  is  almost  two 
magnitude  units.  The  signals  have  the  characteristic  small  upswing,  larger 
downswing  and  largest  second  upswing  observed  in  short  period  explosion  P 
waves  at  teleseismic  distances.  The  usual  nomenclature  for  them  is  the  "a", 
"b"  and  "c"  swings  respectively.  The  "a"  swing  is  very  emergent  at  ALQ  A 
feature  which  could  be  described  as  a  splitting  of  the  "c"  swing  is  indicated 
by  arrows  in  the  figures.  It  is  more  apparent  in  the  higher  frequency  small 
events,  and,  as  is  also  indicated  in  the  figures,  we  believe  it  is  caused  oy 
pP  energy.  One  reason  we  have  for  making  this  inference  is  illustrated  ov 
Figure  3  where  we  show  some  of  the  short  and  long  period  teleseismic  P  waves 
from  the  very  deep  event  CANNIKIN.  As  indicated  by  the  arrows  some  of  the 
data  exhibits  the  same  split  "c"  swing  as  observed  in  the  Pn  data  though  it  is 
washed  out  in  some  cases.  A  synthetic  seismogram  analysis  shows  that  this 
feature  is  caused  by  pP  because  it  appears  when  the  pP  ray  is  added  in 
(Burdick  et  al .  ,  ]dR4).  Synthetics  are  shown  at  the  bottom  where  it  is  also 
illustrated  that  the  variation  in  the  clarity  of  the  feature  could  be  due  to  a 
reasonable  level  of  variation  in  the  level  of  attenuation.  It  could  also  be 
associated  with  azimuthal  variation  in  pP  strength. 

The  implicit  suggestion  from  these  three  figures  that  there  is  a  close 
correlation  between  teleseismic  short  period  P  and  regional  Pn  waveforms  is  a 
significant  one.  It  implies  that  studies  of  the  Pn  waveform  might  be  as 
fruitful  as  studies  of  teleseismic  short  period  P  have  been  for  both 
discrimination  and  yield  estimation  purposes.  The  fundamental  teleseismic 
discriminant,  -  Ms,  relies  in  part  on  short  period  P.  A  related 
discriminant  might  be  built  around  the  first  few  swings  of  Pn  since  that  seems 
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Figure  1.  Observed  Pn  signals  from  Yucca  events  at  DWWSSN  station  ALQ.  Th< 
split  "c"  swing  is  interpreted  as  the  effect  of  pPn.  As  events  get  larger, 
the  increase  in  the  duration  of  the  source  function  washes  out  the  pP  effect. 
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Figure  2.  Observed  P0  signals  from  Yucca  events  at  DWSSN  station  JAS.  There 
Is  a  correspondence  In  the  evolution  of  the  JORNADA  and  ALIGOTE  signals  In 
Figures  1  and  2. 


Figure  3.  The  teleseismic  pP  arrival  can  be  observed  as  a  distinctive  double 
peak  in  the  "c"  swing  of  some  short  period  records  from  CANNIKIN.  This  is  a 
comparison  ^of  observed  short  and  long  period  P  waves  with  synthetics  for  a 
range  of  t  values.  As  t*  increases,  the  interference  in  the  second  upswing 
produced  by  pP  becomes  less  apparent.  The  synthetics  were  computed  with  an 
artificially  delayed  pP. 
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to  be  a  stable  quantity.  There  are  still  many  outstanding  problems  in  the 
area  of  yield  estimation  and  the  understanding  of  explosion  source  physics. 
The  most  successful  previous  seismic  studies  of  these  problems  have  involved 
analysis  of  either  near  field  or  teleseismic  data.  There  is  difficulty  in 
connecting  these  results  because  the  teleseismic  P  waves  are  attenuated  at  an 
unknown  and  frequency  dependent  rate  as  they  travel  through  the  low  velocity 
zone.  The  Pn  waveform  is  not  similarly  affected.  It  can  thus  be  used  to  test 
scaling  laws,  search  for  the  effects  of  spall,  test  models  for  Q,  and  a  number 
of  other  significant  studies.  It  is  interesting  to  compare  the  signals  from 
JORNADA  and  ALIGOTE,  the  largest  and  smallest  events  in  Figures  1  and  2. 
There  is  a  clear  evolution  in  frequency  content  between  the  signals,  and  since 
this  evolution  is  so  similar  at  the  two  stations,  it  must  be  source  and  not 
path  related. 

To  ensure  that  this  splitting  of  the  "c"  swing  is  not  unique  to  the  Yucca 
test  site  or  the  JAS  and  ALQ  stations,  we  also  examined  Pn  data  from  the 
Pahute  test  site  and  from  the  Yucca  site  at  other  stations.  Figures  U  and  5 
show  that  the  Pn  waveform  at  ALQ  and  JAS  for  Pahute  Mesa  events  is  very 
similar  to  that  observed  for  the  Yucca  events.  Another  source  of  easily 
obtainable  data  exists  in  the  form  of  the  analog  WWSSN  film  chips.  In 
reviewing  this  data,  we  found  that  usually  the  gain  was  too  high  to  permit 
inspection  of  the  regional  Pn  waveform  for  NTS  explosions.  However,  whenever 
there  was  an  exceptionally  clear  trace,  some  indication  of  a  split  "cM  swing 
could  be  observed.  Figure  6  shows  some  hand  digitized  traces  from  the  analog 
station  BKS  for  Yucca  events.  The  pP  feature  is  indicated  by  arrows.  Figure 
7  shows  similar  data  for  GSC.  This  data  clearly  shows  that  the  beginning  of 
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PAHUTE  to  ALO  or  AN  MO 
A  =  922  km 
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Figure  4.  Observed  Pn  signals  from  Pahute  events  at  DWWSSN  station  ALQ 
spilt  "c"  swing  Is  Indicated  by  an  arrow. 
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Figure  5.  Observed  Pn  signals  from  Pahute  events  at  DWUSSN  station  JAS. 
split  "c"  swing  is  indicated  by  an  arrow. 


Figure  6.  Observed  Pn  signals  from  Yucca  events  at  WWSSN  analog  station  BKS 
The  split  "c"  swing  is  indicated  by  an  arrow. 


Figure  7. 
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Observed  Pn  signals  from  Yucca  events  at  WWSSN 
"c"  swing  is  Indicated  by  an  arrow. 
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the  P„  waveform  has  a  very  stable  character.  If  it  can  be  shown  that  few 
earthquake  Pn  phases  have  this  character,  we  have  the  basis  of  a  short  period 
regional  discriminant. 

The  similarities  between  teleseismic  short  period  P  waveforms  (Figure  3) 
and  Pn  waveforms  would  indicate  that  there  is  some  similarity  in  the  nature  of 
the  wave  propagation.  If  this  is  true,  then  high  frequency  Pn  should  probably 
not  be  considered  to  be  a  true  head  wave,  but  rather  a  turning  ray  in  the 
smooth  gradient  of  the  lid.  To  investigate  the  character  of  regional  Pn  in 
the  western  U.  S.  in  more  detail,  we  measured  its  speed  at  each  of  the  four 
stations  discussed  above.  The  results  are  summarized  in  Figure  8. 
Time-distance  plots  are  shown  for  each  station  along  with  a  best  fitting  line 
and  the  apparent  velocity.  Wave  speed  does  not  increase  monotonically  with 
distance  because  the  BKS  speed  is  slower  than  the  JAS.  This  only  indicates, 
however,  that  velocity  varies  laterally  in  the  lid.  Furthermore,  there  is 
enough  of  a  gradient  in  velocity  (6.84  to  8.23  km/sec  between  2  and  8  degrees) 
to  support  a  strong  turning  ray.  A  similar  increase  in  apparent  velocity  is 
associated  with  teleseismic  rays  spanning  the  ranges  between  30  and  45°.  It 
is  therefore  not  unreasonable  to  suggest  that  Pn  is  a  turning  ray  from  the  lid 
between  2  and  10°,  at  least  in  the  western  U.  S. 

For  our  initial  calculations,  we  assume  that  the  turning  ray  model  for 
Pn  is  appropriate.  We  utilize  the  same  methodology  and  assumptions  used  for 
computing  teleseismic  short  period  P.  The  response  of  the  earth  is  assumed  to 
be  a  delta  function  with  an  amplitude  controlled  by  geometric  spreading  with 
the  exception  of  the  effects  of  near  source  structure.  To  simulate  this 
effect,  we  use  a  plane  wave  calculation  performed  using  propagator  matrices. 
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Figure  8.  Determinations  of  wave  speed  of  Pn  at  regional  stations  around  NTS. 
Such  information  would  be  routinely  available  from  an  intelligent  array  and 
could  be  used  in  a  discrimination  algorithm. 


We  assumed  a  near  source  crustal  structure  based  on  modeling  of  near  field 
data  (Hartzell  et  al.,  1983),  and  used  a  Pahute  scaled  Mueller -Murphy  source. 
An  apparent  velocity  of  8.3  km/sec  was  used  which  is  appropriate  for  ALQ.  The 
synthetics  for  four  different  source  depths  are  compared  to  observed  Pahute 
Mesa  Pn  phases  in  Figure  9.  The  observed  feature  in  the  "c"  swing  of  the  data 
is  notably  later  in  the  observations  than  the  synthetics,  but  the  match  of  the 
bottom  traces  is  actually  quite  good.  The  observation  that  actual  pP  arrival 
time  is  later  than  one  would  expect  from  elastic  predictions  has  been  made 
commonly  (Burdick  et  al.,  1984a, b)  and  is  generally  interpreted  to  mean  that 
the  arrival  is  a  combination  of  elastic  and  inelastic  effects  over  the  source. 
Although  we  refer  to  the  feature  as  pP  here  we  acknowledge  that  it  is  probably 
something  more  complex.  The  key  point,  however,  is  that  it  is  a  depth  related 
phenomenon  and  something  thus  could  cause  the  basis  for  a  discriminant.  A 
correction  for  attenuation  was  initially  included  in  the  calculations. 
However,  we  found  that  the  periodicity  of  the  data  was  best  matched  for 
vanishingly  small  t*.  Q  for  Pn  is  apparently  very  high  in  the  western  U.  S. 

As  was  noted  earlier,  the  "a"  swing  at  ALQ  is  unusually  broad  and  low  in 
relative  amplitude.  We  indicate  the  beginning  of  the  "a"  swing  with  an  arrow 
in  Figure  9.  If  Pn  behaves  as  a  true  head  wave,  it  should  begin  more 
gradually  than  if  it  were  a  turning  ray.  Thus,  we  may  be  observing  some 
indication  that  this  is  the  case  at  ALQ.  At  high  frequency,  an  idealized  head 
wave  can  be  well  approximated  as  the  integral  of  the  turning  ray.  We  compare 
integrated  synthetics  to  the  ALQ  observations  in  Figure  10.  The  smoothing 
effect  can  be  seen  in  the  differences  between  Figures  9  and  10,  but  the 
extreme  smoothing  of  just  the  "a"  swing  is  not  matched.  The  pP  feature  begins 
to  emerge,  but  not  as  strongly  as  for  the  turning  ray  case.  It  presumably 
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Figure  9.  A  comparison  of  synthetic  to  observed  ALQ  records  from  NTS 
explosions.  We  assume  that  PD  is  a  simple  turning  ray  in  the  lid  and  perform 
a  plane  wave  calculation  using  a  structure  for  Pahute  Mesa.  The  effective  pP 
is  obviously  late. 


would  for  a  deep  enough  depth,  but  the  synthetics  are  already  too  long  in 
period  to  match  the  data  for  the  deepest  depth  shown.  We  believe  that  the 
representation  of  high  frequency  Pn  as  a  turning  ray  is  more  appropriate.  The 
unusual  behavior  at  ALQ  is  more  likely  related  to  lateral  variations  in 
structure  or  to  the  frequency  dependence  of  Q.  In  the  latter  case,  we  could 
be  seeing  some  form  of  dispersion  where  the  highest  frequency  energy  arrives 
measurably  ahead  of  the  rest  causing  a  precursor. 

Averaging  :  We  have  shown  in  previous  figures  that  there  can  be  little 
doubt  that  the  splitting  of  the  "c"  swing  by  pP  is  a  fairly  general 
occurrence.  There  are,  however,  variations  in  the  waveshape  at  each  station. 
Some  appear  to  be  systematically  related  to  event  size  and  others  to  random 
scattering.  In  this  section,  we  utilize  each  station  as  an  array  to  stack  for 
the  average  properties  of  the  wave  shape.  The  results  of  the  averaging  for 
the  four  stations  are  shown  in  Figures  11  and  12.  The  traces  have  been 
aligned  to  preserve  the  splitting  of  the  "c"  swing.  The  average  Pn  phases  are 
shown  together  in  Figure  13.  The  pulse  shapes  are  most  similar  at  GSC  and 
BKS.  JAS  seems  to  have  a  much  more  deeply  split  "c"  swing,  but  it  otherwise 
correlates  well  with  BKS.  The  pP  feature  appears  to  arrive  earlier  in  the 
record  from  the  most  distant  station,  ALQ,  as  one  should  expect  to  see.  It  is 
interesting  to  note  that  the  average  Pn  waveform  does  not  change  at  a  fixed 
station  for  the  two  different  test  sites.  The  average  Pn  waveform  at  JAS  and 
ALQ  is  the  same  for  Yucca  and  Pahute  even  though  the  test  sites  are  more  than 
50  km  apart.  This  implies  that  the  average  Pn  waveforms  at  the  stations  are 
stable  and  not  strongly  affected  by  near  source  structure.  This  is  an 
important  fact  because  the  ability  to  discriminate  based  on  this  pP  phenomenon 
will  depend  strongly  on  how  well  the  waveform  for  a  new  source-site 
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Figure  11.  The  average  Pn  waveshape  at  CSC  and  JAS  is  relatively  consistent 
and  displays  the  split  "c"  swing.  The  sum  of  the  individual  traces  is  shown 
at  the  bottom. 
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Figure  12.  The  average  Pn  waveshape  at  BKS  and  ALQ  Is  relatively  consistent 
and  displays  the  split  "c"  swing.  The  sum  of  the  individual  traces  Is  shown 
at  the  bottom. 
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Figure  13.  A  summary  of  the  average  Pn's  regional  from  stations  around  NTS 
The  average  Pn  does  not  vary  at  a  given  station  between  Pahute  and  Yucca. 


combination  can  be  predicted.  The  conclusion  of  this  exercise  is  that  the  pP 
interference  is  a  pervasive  feature,  and  insofar  as  it  really  is  caused  by  the 
free  surface,  one  should  expect  to  see  it  for  any  contained  test.  The  next 
question  which  must  be  answered  is  how  often  similar  features  appear  in 
earthquake  waveforms. 

Discrimination  using  Pj,:  To  address  this  question,  we  assembled  a  data  base 
of  31  observations  of  western  United  States  earthquakes  from  stations  ALQ  and 
JAS.  The  earthquakes  were  selected  to  be  as  close  to  NTS  as  possible,  but  to 
obtain  enough  samples  we  needed  to  accept  events  from  as  far  away  as  southern 
California.  We  estimated  the  Pn  arrival  time  as  accurately  as  we  could  and 
windowed  out  the  first  three  seconds  of  the  arrival  Representative  samples 
which  were  selected  at  random  are  shown  in  Figures  14  and  15.  We  compared 
these  signals  to  the  average  explosion  Pn  waveforms  shown  in  Figure  12  using 
the  analytic  norm  introduced  by  Burdick  and  Mellman  (1976).  This  norm 
increases  uniquely  to  a  value  of  one  if  and  only  if  two  signals  are  identical. 
It  can  be  shown  to  be  closely  related  to  the  L2  norm  for  two  traces  if  they 
are  aligned  at  the  time  of  maximum  correlation  (Burdick,  1985). 

The  results  for  Pahute  events  at  JAS  are  shown  in  Figure  16.  The 
separation  of  the  explosion  and  earthquake  populations  is  complete  indicating 
that  the  beginning  of  an  earthquake  signal  is  systematically  different  from 
the  beginning  of  an  explosion  signal.  The  results  for  Pahute  events  at  ALQ 
are  shown  in  Figure  17.  Again  the  discrimination  between  events  is  quite 
good.  Unfortunately,  most  events  from  Pahute  Mesa  are  relatively  large.  The 
discrimination  shown  in  Figures  16  and  17  is  not  in  the  small  (nij,  <  4.0) 
magnitude  range  where  discrimination  is  difficult  using  standard  techniques. 
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RANDOMLY  SELECTED 
EARTHQUAKE  Pn  ARRIVALS  AT  JAS 
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Figure  15.  A  data  base  of  digital  PB  arrivals  from  western  U.  S.  earthquakes 
has  been  assembled  from  station  JAS.  These  are  the  first  6  arrivals  chosen. 
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Figure  16.  Discrimination  of  Pahute  explosions  from  earthquakes  using  the 
correlation  with  the  average  Pn  waveform  at  JAS.  The  explosions  are  stars  and 
the  earthquakes  squares. 
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Figure  17.  Discrimination  of  Pahute  explosions  from  earthquakes  using  the 
correlation  with  the  average  Pn  waveform  at  A1Q.  The  explosions  are  stars  and 
the  earthquakes  squares. 
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Smaller  events  are  generally  located  at  Yucca  Flat.  The  results  for 
Yucca  events  at  JAS  are  shown  in  Figure  18.  The  stars  were  generated  using 
the  Pn  waveform  shown  in  Figure  11.  The  separation  is  again  good  above 
magnitude  5  but  definitely  begins  to  deteriorate  as  magnitude  decreases  to  3. 


1 
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3), 

and 

its 

1 

correlation 

with  the  average 
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Nonetheless,  the 

Pn  waveform 

does  appear  to 

have 
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discrimination  power  between  magnitude  3  and  U  where  new  discriminants  are 
needed. 


Part  of  the  deterioration  of  the  discriminant  at  small  magnitude  is  due 
to  the  fact  that  the  average  waveform  used  for  the  correlations  came  from 
larger  events.  We  averaged  the  waveforms  from  6  small  explosions  and 
correlated  them  with  the  earthquakes  to  produce  the  crosses  in  Figure  18 .  The 
event  PINEAU  is  clearly  separated  from  the  earthquake  population.  The  large 
explosions  are  not  nearly  as  well  correlated  with  this  estimate  of  average  Pn 
(i.e.  the  crosses  are  lower  than  the  stars  at  large  magnitude)  indicating  that 
there  is  a  size  effect  which  needs  to  be  taken  into  account  in  using  this 
discriminant . 


Figure  19  shows  the  results  for  Yucca  events  at  ALQ.  Again  the 
populations  are  separated,  but  discrimination  power  decreases  with  magnitude. 
We  have  examined  some  of  the  worst  cases  and  find  that  this  is  primarily  a 
signal  to  noise  effect.  It  is  important  to  note  that  if  we  were  working  with 
array  data  like  that  from  NORESS,  we  could  process  the  signals  to  improve 
signal  to  noise.  Moreover,  it  is  well  known  that  Pn  amplitudes  are  generally 
much  higher  in  stable  continental  areas  than  the  western  U.  S.,  so  we  can 
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Figure  18.  Discrimination  of  Yucca  explosions  from  earthquakes  using  the 
correlation  with  the  average  Pn  waveform  at  JAS.  The  explosions  are  stars 
and  crosses  and  the  earthquakes  are  squares  and.  The  stars  represent 
correlation  with  the  average  moderate  event  waveform  and  the  crosses  with  the 
average  small  event. 
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Figure  19.  Discrimination  of  Yucca  explosions  from  earthquakes  using  the 
correlation  with  the  average  Pn  waveform  at  ALQ.  The  explosions  are  stars  and 
the  earthquakes  are  squares. 


anticipate  good  signal  to  noise  ratio  down  to  magnitude  3  and  less. 

As  a  final  exercise,  we  test  the  observation  made  earlier  in  Figure  13 
that  there  does  not  seem  to  be  much  variation  in  average  Pn  waveform  between 
the  two  test  sites  at  a  fixed  station  using  the  quantitative  norm.  In  Figures 
20  and  21,  we  show  the  correlation  between  all  explosion  waveforms  and  the 
average  Yucca  Pn  waveform  for  ALQ  and  JAS  respectively.  There  is  effectively 
no  separation  of  the  two  populations  indicating  that  the  average  wave  shape  is 
not  strongly  dependent  on  near  source  structure.  This  is  a  result  that  is 
surprising  at  such  high  frequencies  but  welcome  because  it  indicates  that  we 
may  be  able  to  develop  a  stable  discriminant. 

Such  a  discriminant  could  be  based  on  a  wide  variety  of  underlying 
strategies  and  waveform  comparators  besides  the  one  we  have  been  using  so  far. 
The  important  result  we  have  shown  is  just  that  the  onset  of  Pn  for  explosions 
is  stable  and  unique.  One  strategy  for  a  discrimination  procedure  would  be  a 
purely  empirical  one.  As  large  a  data  base  of  average  Pn  waveforms  as 
possible  could  be  assembled  from  U.  S.,  Soviet  and  French  test  sites.  This 
data  base  could  then  be  compared  continuously  to  the  digital  Pn  signals  as 
they  stream  into  the  global  array.  With  modern  processors,  this  operation 
could  be  done  very  rapidly  for  a  large  background  data  base  of  average  Pn 
phases.  The  decision  could  continuously  be  made  whether  Pn  signals  from  new 
events  look  like  Pn  phases  observed  from  previous  explosions.  An  alternate 
method  could  involve  development  of  a  synthetic  predictor  of  Pn  waveforms 
based  on  source  to  station  distance,  event  magnitude,  observed  Pn  wavespeed 
and  other  source  parameters.  A  third  approach  could  involve  the  new  computer 
based  recognition  techniques  highly  tuned  to  the  Pn  waveform. 
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Figure  20.  Correlation  of  Pn  waveforms  from  either  Pahute  or  Yucca  events 
with  the  average  explosion  Pn  waveform  at  ALQ. 
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Figure  21.  Correlation  of  PB  waveforms  from  either  Pahute  or  Yucca  events 
with  the  average  explosion  Pn  waveform  at  JAS . 


The  Crustal  Resonance  Phase  Discriminant 


The  chief  drawback  to  using  regional  Pn  as  the  basis  for  a  discriminant  | 

is  that  it  does  not  have  as  high  an  amplitude  as  other  regional  phases.  Its 
advantage  is  that  it  is  the  first  arrival  on  the  record,  and  the  idea  of  a 
waveform  for  it  is  meaningful  as  was  demonstrated  in  the  preceding  section.  | 

In  this  section,  we  discuss  an  attempt  to  base  a  discriminant  on  some  of  the 
higher  amplitude  arrivals.  The  phase  Pg  was  selected  for  study  in  this,  the 
developmental  phase  of  the  project,  where  data  from  the  western  U.  S.  are  | 

being  analyzed.  The  reason  for  making  this  selection  is  illustrated  in  Figure 
22.  A  suite  of  digital  regional  waveforms  from  explosions  and  earthquakes  are 
displayed  in  a  seismic  section.  In  terms  of  signal  to  noise  ratio,  the  most  | 

effective  discriminants  should  be  based  either  on  Pg  or  the  Sg  -  Lg  complex. 

We  indicate  the  predicted  moveout  of  PmP,  Sn  and  SmS  for  an  appropriate 

crustal  model  to  illustrate  two  points.  The  first  is  that  the  Pn  -  Sn  I 

discriminant  which  was  developed  using  NORESS  data  by  Dysart  and  Pulli  (1988) 

can  not  be  tested  in  this  tectonic  environment.  There  is  effectively  no  Sn  in 

this  frequency  band.  The  second  point  is  that  the  Pg  and  Sg  wave  complexes 

begin  at  PmP  and  SmS  times  respectively  although  they  have  a  duration  many 

times  longer  than  could  be  explained  by  just  these  phases.  The  reason  for  the 

long  duration  of  at  least  Pg  and  possibly  Sg  as  well  was  discussed  by  Burdick 

and  Helmberger  (1988)  who  show  that  the  complex  is  made  up  of  a  series  of 

pulses  called  crustal  resonance  phases.  They  are  built  up  by  groups  of  rays 

reverberating  in  the  crust.  This  sequence  of  pulses  and  their  coda  account 

for  the  character  of  Pg  Similar  rays  in  the  shear  mode  may  account  for  Sg. 

The  advantage  of  studies  made  in  the  time  domain  is  that  deterministic 
concepts  seem  to  apply  much  better  at  times  before  coda  builds  up.  It  seems 
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Figure  22.  A  short  period  seismic  section  of  explosions  and  earthquakes  for 
the  western  U.  S,  The  Pg  and  S  arrivals  commence  with  PmP  and  SmP.  There  is 
no  3n  arrival  in  this  frequency  band. 
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logical  therefore  to  analyze  the  beginning  of  Pg  to  find  any  deterministic 
character  on  which  a  discriminant  could  be  based  since  it  is  the  first  of  the 
high  amplitude  regional  arrivals  and  it  generally  overwhelms  Pn. 

There  have  been  several  deterministic  modeling  studies  of  the  onset  of  Pg 
although  they  have  not  widely  been  recognized  for  what  they  are.  This  is 
because  they  have  been  performed  on  data  in  the  long  period  band  and  are 
generally  termed  Pnl  studies.  Though  the  emphasis  in  these  studies  was  on  the 
long  period  energy,  the  forward  synthetics  often  match  the  short  period 
details  of  observed  seismograms  remarkably  well.  To  illustrate  the 
correspondence  between  Pg  and  the  arrivals  in  the  long  period  band,  we  show  a 
section  of  long  period  records  from  NTS  explosions  in  Figure  23  over  a  range 
comparable  to  that  shown  in  Figure  22.  The  long  period  energy  at  the 
beginning  of  the  records  is  made  up  of  a  complex  of  head  waves  and  the  shift 
to  higher  frequency  is  related  to  a  complex  of  reflected  rays.  The 
correspondence  of  this  reflected  energy  and  Pg  is  obvious.  The  data  in  Figure 
23  were  the  data  used  in  a  well  known  study  by  Wallace  et  al.  (1983)  to 
deterministically  characterize  the  effects  of  tectonic  release.  We  have  also 
deterministically  modeled  it  along  with  comparable  western  U.  S.  earthquake 
data  as  described  in  a  previous  report  on  this  study  (Burdick  and  Helmberger, 
1988)  . 

To  further  this  phase  of  the  work  and  to  ensure  that  the  results  from  the 
western  U.  S.  will  be  applicable  in  other  tectonic  environments,  we  have 
examined  Pnl  waves  from  two  eastern  U.  S.  earthquakes  recorded  at  eastern  U.S. 
and  Canadian  stations.  These  were  the  9/4/63  Baffin  Bay  earthquake  and  the 
1/9/82  New  Brunswick  event.  The  results  are  shown  in  Figures  24  and  25.  As 
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Figure  23.  A  long  period  seismic  section  of  explosion  records  from  the 
western  U.  S.  The  high  frequency  energy  associated  with  the  Pg  phase  is  clear 
in  both  Figure  22  and  23  and  commences  with  the  PmP  arrival. 


usual,  the  long  periods  are  fit  extremely  well,  but  so  are  many  of  the  details 
at  short  periods.  The  signal  appears  deterministic  down  to  periods  perhaps  as 
short  as  two  seconds.  We  are  not  suggesting  that  this  deterministic  behavior 
holds  up  at  high  frequency,  but  we  do  suggest  that  Pnl  type  calculations  might 
be  used  as  a  guide  in  identifying  when  information  that  is  rich  in 
discrimination  potential  is  arriving  in  the  record.  We  also  conclude  that  our 
western  U.  S.  studies  of  crustal  resonance  phases  will  be  relevant  to  data 
from  other  tectonic  environs. 

Crustal  Resonance  Phases.:  It  would  be  useful  at  this  juncture  to  review  some 
of  the  details  associated  with  computation  and  analysis  of  crustal  resonance 
phases.  Just  as  in  Pnl  calculations,  we  assume  a  layer  over  a  half  space 
crustal  model.  A  generalized  ray  representation  of  the  response  of  the  model 
is  particularly  well  suited  for  determining  exactly  how  arrivals  such  as 
crustal  resonance  phases  are  being  excited.  The  Green's  function  can  be  built 
progressively  by  adding  in  more  rays  until  the  feature  of  interest  emerges. 
Figure  26  shows  the  development  of  the  Green's  function  for  the  Helmberger  and 
Engen  (1980)  western  U.  S.  crust  model  at  a  distance  of  1000  km.  It  is  being 
progressively  constructed  through  sequential  ray  addition.  Two  of  the  four 
fundamental  Green's  functions  are  shown;  vertical  strike  slip  on  the  left  and 
vertical  dip  slip  on  the  right.  These  two  fundamental  Green's  functions 
produce  the  most  variable  regional  P  waves  because  they  have  the  strongest 
variation  in  vertical  radiation  pattern  as  shown  at  the  top.  The  response  to 
just  the  direct  P  and  S  waves  is  shown  first.  Then  the  effect  of  adding  in 
those  rays  with  just  one  reverberation  in  the  crust  along  with  their 
associated  depth  phases  (PmP,  pPmP,  SmP,  sSmP  etc.)is  shown.  The  third 
response  is  for  two  reverberations  in  the  crust  with  depth  phases,  the  fourth 
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Figure  24.  Observed  and  synthetic  Pnl  waveforms  for  the  September  4 
Baffin  Bay  earthquake.  Peak  to  peak  amplitudes  in  cm  x  10‘3  are  shown. 
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Figure  26.  Green's  functions  for  a  layer  over  a  half  space  computed  using 
generalized  rays.  Successive  rows  are  computed  by  adding  in  more  crustal 
reverberations  as  indicated.  The  crustal  resonance  phases  are  high  frequency 
arrivals  commencing  about  half  way  through  the  trace. 
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for  four  and  so  on.  The  number  of  additional  rays  in  each  iteration  is 
indicated.  The  number  of  rays  grows  so  rapidly  with  each  iteration  that 
considering  any  more  than  5  reverberations  would  be  impractical.  The  crustal 
resonance  phases  are  the  periodic  high  frequency  pulses  commencing  about  half 
way  through  the  trace.  Successive  resonance  phases  appear  as  successive 
reverberations  of  rays  are  added  in.  The  long  period  oscillations  that  build 
up  have  classically  been  called  PL. 

Since  it  is  clear  that  these  reverberation  phases  exist,  at  least  at 
moderate  periods,  it  is  natural  to  explore  the  possibility  that  they  could  be 
useful  in  discrimination.  The  use  of  ratios  of  phases  for  discrimination  such 
as  ratios  of  Pg  to  Lg,  suggest  a  similar  study  of  the  ratio  of  the  first  to 
the  second  crustal  resonance?  In  further  work,  we  examined  the  composition  of 
the  first  two  resonance  phases  in  terms  of  generalized  rays  in  more  detail. 
For  an  explosion  the  first  resonance  is  almost  completely  made  up  of  the  three 
rays  P,  PmP  and  pPmP.  Its  high  frequency  is  caused  by  the  strong  interference 
of  PmP  with  the  depth  phase.  The  direct  ray,  P,  contributes  primarily  long 
period  energy.  Interestingly,  the  second  resonance  phase  is  not  dominated  by 
rays  with  a  complete  extra  reverberation  in  the  crust,  but  by  rays  which  have 
at  least  one  converted  leg.  For  an  explosion,  the  rays  are  PmS,  pSraS,  PmPSmS, 
PmPSmP,  and  PmSPmP.  Earthquakes  would  also  contribute  energy  to  the  second 
resonance  through  SmS,  sPmP  and  SmSPmP.  However,  it  is  important  to  bear  in 
mind  that  this  is  high  frequency  regional  data,  and  only  stable  properties  of 
the  data  are  important.  As  we  will  show  in  a  later  section,  the  first  crustal 
resonance  appears  to  be  an  efficient  energy  channel  for  P  and  the  second 


resonance  for  S.  The  best  method  for  discriminating  using  Pg  may  be  to  look 
at  the  relative  amount  of  scattered  SV  energy  between  the  first  and  second 
resonance  phases.  There  should  be  more  scattered  S  in  earthquakes. 

The  basic  philosophy  behind  many  discrimination  studies  has  been  to 
search  for  stronger  S  excitation  than  P  for  earthquakes  in  one  pair  of 
regional  phases  or  another.  The  best  known  studies  are  probably  those 
involving  the  ratio  of  Pg  to  Lg  or  Pn  to  Sn  (Blandford,  1981;  Dysart  and  Pulli, 
1988).  To  those  familiar  with  that  work,  it  may  seem  strange  to  consider 
looking  at  the  relative  P  and  S  excitation  within  just  Pg.  However,  there  are 
strong  channels  for  energy  that  leaves  the  source  as  S  into  different  crustal 
resonances  and  they  differ  from  resonance  to  resonance.  It  is  reasonable  to 
expect  a  change  in  the  amount  of  S  energy  between  the  first  and  second 
resonance  and  that  the  ratio  of  this  energy  will  differ  between  explosions  and 
earthquakes . 

There  is  no  doubt  that  crustal  resonance  phases  exist  in  an  idealized 
case.  However,  it  remains  to  be  established  to  what  extent  they  exist  and  how 
coherent  they  are  in  a  real  crust  at  high  frequency.  To  search  for  these 
phases,  we  used  the  same  digital  data  base  from  ALQ  and  ANMO  that  we  used  in 
the  Pn  study  discussed  in  the  preceding,  but  concentrated  on  the  Pg  portion  of 
the  signal.  Although  there  were  some  indications  of  a  series  of  arrivals 
pumping  energy  into  the  Pg  coda  in  some  of  the  raw  data,  the  effects  of 
scattering  and  other  variations  between  events  seems  to  dominate.  This  has 
long  been  known  to  be  the  greatest  difficulty  in  interpreting  regional  data. 
To  suppress  the  scattering  effects,  we  aligned  the  signals  very  carefully  at 
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Pn  time  and  summed.  The  Yucca  and  Pahute  data  were  treated  separately.  This 
procedure  neglects  the  difference  in  slowness  between  resonance  phases  and  Pn, 
but  our  results  indicate  that  this  is  a  reasonable  approach. 

The  results  are  shown  in  Figures  27  and  28  where  we  display  progressive 
sums  of  records.  Each  subsequent  line  shows  the  effects  of  adding  in  an 
additional  record.  The  final  sums  clearly  indicate  that  there  are  coherent 
resonance  phases  at  frequencies  of  several  hz.  There  are  10  traces  summed  in 
the  stack  of  Yucca  events  in  Figure  27.  The  first  two  crustal  resonance 
phases  develop  almost  immediately  indicating  that  they  are  very  coherent.  A 
third  appears  by  the  time  4  traces  are  added  in  and  the  sum  has  stabilized  by 
the  time  8  are  added  in.  Actually,  this  is  a  relatively  rapid  rate  of 
stabilization  indicating  that  the  resonance  phases  are  very  coherent  for  this 
path.  In  the  Pahute  record  stack,  the  first  and  third  resonance  phases 
stabilize  first,  while  the  second  does  not  really  emerge  clearly  until  11 
traces  are  added.  At  the  bottom,  we  show  a  synthetic  for  a  simple  layer  over 
a  half  space  model  just  to  show  that  the  resonance  phases  have  approximately 
the  character  we  expect. 

An  important  instrumentation  change  occurred  at  the  ANMO  station  in  about 
1983.  Two  horizontal  short  period  instruments  were  added  which  began  to 
supply  the  necessary  information  for  studying  the  relative  polarization  of 
crustal  resonances.  A  limited  number  of  explosions  have  been  recorded  so  far. 
We  have  performed  some  preliminary  analyses  of  this  data.  By  simply 
inspecting  the  data,  we  have  developed  the  impression  that  the  first  resonance 
is  an  efficient  compressional  energy  channel  and  the  second  is  much  more 
favorable  for  shear.  Figure  29  shows  a  good  example  of  the  effect  we  are 
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Figure  27.  The  crustal  resonance  phases  for  Yucca  events  are  separated  from 
the  noise  by  stacking.  The  synthetic  at  the  bottom  merely  indicates  that  the 
phases  have  approximately  the  right  timing  and  amplitude. 
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Figure  28.  The  crustal  resonance  phases  for  Pahute  events  are  separated  from 
the  noise  by  stacking.  The  synthetic  at  the  bottom  merely  indicates  that  the 
ph.i-es  have  approximately  the  right  timing  and  amplitude. 
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seeing.  In  this  particular  instance,  the  crustal  resonances  are  visible 
without  stacking.  They  appear  to  be  very  similar  on  the  radial  and  vertical. 
The  second  resonance  is  heavily  favored  on  the  tangential,  however.  This 
could  not  be  true  SH  energy  in  the  theoretical  sense,  but  recalling  that  this 
is  high  frequency  energy  traveling  in  the  real  earth  it  could  be  interpreted 
as  scattered  S  traveling  in  one  or  several  of  the  favored  shear  modes. 

The  majority  of  the  events  for  which  we  have  data  so  far  are  from  Pahute 
since  those  events  are  larger  and  trigger  the  instrument  more  consistently. 
We  have  rotated  and  stacked  the  traces  of  10  events.  The  results  are  shown  in 
Figure  30.  On  average,  the  second  resonance  is  suppressed  on  the  radial 
relative  to  either  the  vertical  or  tangential.  The  first  is  equally  strong  on 
all  three.  This  is  not  inconsistent  with  the  example  shown  in  the  previous 
figure.  On  average,  there  could  be  a  component  of  off  azimuth  P  on  the 
tangential  to  account  for  the  first  arrival  and  then  a  consistent  component  of 
scattered  S  to  account  for  the  second.  The  suppression  of  the  second 
resonance  on  the  radial  component  is  interesting  because  this  is  what  is 
predicted  by  the  synthetics  shown  on  the  bottom.  The  SV  energy  must  approach 
the  station  at  low  apparent  velocity  and  be  concentrated  on  the  vertical.  The 
key  point  is  that  there  is  a  clear  evolution  of  particle  motion  between  the 
first  and  second  resonances.  The  questions  ahead  are  can  it  be  reliably 
measured  for  single  events  and  can  it  be  used  to  discriminate? 

Another  important  type  of  information  that  can  be  resolved  through  the 
stacking  procedure  illustrated  by  Figure  27  regards  the  character  of  the 
signal  generated  noise  in  this  type  of  data.  The  traces  are  all  normalized  to 
unit  peak  before  stacking.  The  peak  amplitudes  of  the  sums  after  stacking  are 
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Figure  30.  Average  three  component  Pg  waveforms  for  Pahute  events  at  ANMO  are 
determined  by  stacking.  Note  the  evolution  of  particle  motion  between  the 
first  and  second  crustal  resonance.  A  similar  evolution  is  predicted  in  the 
synthetics . 


indicated.  The  average  peak  Tangential/Radial  prior  to  stacking  was  .515. 
After  stacking  it  was  only  reduced  to  .351.  There  is  a  significant  component 
of  the  energy  in  the  P  wave  train  that  is  both  coherent  and  tangential.  This 
type  of  noise  data  could  be  used  in  predicting  the  expected  performance  of 
discriminants  that  should  work  in  theory. 


DISCUSSION 


There  is  little  need  to  discuss  the  Pn  discriminant  further.  It  is  simple 
in  concept  and  design  which  makes  the  likelihood  that  it  will  work  relatively 
high.  The  basic  idea  behind  it  is  that  at  short  times,  before  the  depth 
phases  from  an  earthquake  could  possibly  arrive,  the  explosion  waveform 
differs  systematically  from  the  earthquake  because  its  depth  phases  have 
arrived. 

The  Pg  discriminant  is  in  a  much  earlier  phase  of  development  than  the  Pr , 
but  the  stacking  experiments  show  that  there  is  a  strong  basis  for  further 
work.  Deterministic  calculations  in  highly  idealized  crust  models  have  played 
an  important  role  in  the  analysis  of  the  discriminant  presented  here.  It  is 
important  to  note  that  it  is  not  intended  that  such  calculations  will  play  a 
direct  role  in  the  actual  application  of  the  discriminant.  They  are  now  only 
being  used  as  a  tool  for  interpretation  of  the  stacked  data  to  help  guide  in 
choosing  which  features  of  regional  arrivals  are  most  likely  to  differ  between 
earthquakes  and  explosions. 

It  would  seem  that  most  researchers  who  have  dealt  with  high  frequency 
regional  data  have  rapidly  abandoned  a  forward  modeling  approach  to  data 
interpretation  and  taken  a  purely  empirical  approach.  There  are  many  good 
reasons  for  this.  The  basic  principles  which  apply  at  periods  longer  than  a 
second  clearly  break  down  in  a  severe  way  at  higher  frequency.  The  duration 
of  regional  phases  is  much  longer  than  the  possible  duration  of  the  source, 
and  coupling  between  the  radial  and  tangential  components  is  very  strong.  In 
this  work,  we  have  shown  through  st<  eking  of  the  Pg  complex  that  on  average  it 
is  made  up  of  a  series  of  coherent  subpulses  which  individually  do  have 
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reasonable  durations.  The  synthetic  seismogram  modeling  allows  interpretation 


of  the  pulses  as  crustal  resonance  phases.  The  continuing  arrival  of  these 
pulses  along  with  a  reasonable  amount  of  coda  account  for  the  overall  duration 
of  the  Pg  complex.  Thus,  though  the  calculations  are  highly  idealized,  they 
have  provided  an  important  guide  to  interpretation.  The  calculations  also 
indicate  that  there  are  strong  channels  for  S  into  the  Pg  phase  and  for  P  into 
the  Sg  -  Lg  complex.  These  channels  may  account  for  why  the  Pg  to  Lg 
discriminant  performs  poorly.  However,  it  would  also  appear  that  the  ratios 
of  the  first  to  second  resonance  phase  may  constitute  the  basis  for  a 
discriminant .  No  forward  synthetics  would  be  involved  in  the  actual 
implementation  of  the  discriminant.  This  would  require  sophisticated  phase 
detectors  and  particle  motion  filters  which  would  measure  average  properties 
of  Pg  with  time . 
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CONCLUSIONS 


The  overall  purpose  of  this  work  is  to  develop  and  test  regional 
discriminants  in  the  western  U.  S.  and  then  to  investigate  how  reliably  they 
can  be  transported  to  other  areas  where  earth  structure  and  tectonic 
conditions  might  be  much  different.  The  Pn  waveform  discriminant  has  been 
developed  in  principle.  Future  efforts  need  to  be  directed  towards  measuring 
additional  average  Pn  waveforms  at  digital  stations  in  the  western  U.  S.  The 
Pg  discriminant  is  in  a  much  more  primitive  state.  Additional  studies  are 
needed  to  fully  characterize  crustal  resonance  phases  with  particular 
attention  being  paid  to  relative  particle  motion.  The  studies  of  Pnl  in  the 
eastern  U.  S.  and  Canadian  shield  indicate  that  these  resonances  are  pervasive 
phenomena,  and  if  a  discriminant  can  be  developed,  it  should  be  highly 
transportable.  We  note,  in  this  regard,  that  Vogfjord  and  Langston  (1988) 
have  already  reported  the  detection  of  what  are  probably  resonance  phases  at 
NORESS.  They  called  them  Pgl  and  Pg2  however. 
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